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ABSTRACT 

Fluctuating  velocity  and  temperature  data  collected  from  instru- 
ments mounted  at  three  levels  on  the  R/V  ACANIA  were  analyzed.   The 
structure  of  atmospheric  turbulence  was  examined  on  the  basis  of 
spectral  distributions  of  wind  and  temperature  fluctuations  over  the 
ocean.   The  results  of  this  study  were  utilized  to  determine  the 
validity  of  present  formulations  and  prediction  techniques  for  the 
momentum  flux  and  the  temperature-structure  parameter  over  water. 

For  the  stability  conditions  encountered,  the  momentum  flux 
calculations  support  present  expressions  for  the  constant  flux  layer. 
Proper  directional  deviation  from  prediction  of  the  variation  of  e 
with  height  was  observed.   When  corrections  for  stability  were  applied, 
proper  adjustment  toward  the  predictions  were  noted.   However,  com- 
plete adjustment  was  not  evident.   This  suggests  other  factors  as 

2 
wind  wave  effects  may  be  important.   C    calculations  were  in  agree- 
ment with  present  formulations  and  prediction  techniques. 
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I.   INTRODUCTION 

Interdisciplinary  investigations  are  being  conducted  at  the  Naval 
Postgraduate  School  to  correlate  optical  propagation  and  atmospheric 
turbulence  data  in  the  marine  boundary  layer.   Such  information  is 
necessary  for  design  of  future  naval  systems.   In  these  investigations, 
personnel  from  the  Meteorology,  Mechanical  Engineering,  Oceanography, 
and  Physics  departments  have  been  obtaining  shipboard  observational 
data  to  evaluate  optical  propagation  in  the  near-ocean  environment. 

The  purpose  of  this  study  is  to  examine  the  structure  of  atmos- 
pheric turbulence  on  the  basis  of  spectral  distributions  of  wind  and 
temperature  fluctuations.   Turbulence  similarity  expressions  have  been 
tested  primarily  over  land.   Only  a  few  over  water  results  have  been 
obtained.   Observational  results  from  this  study  will  provide  a  basis 
on  which  to  evaluate  and  improve  existing  techniques  in  predicting 
over  water  optical  propagation. 

During  January  1974,  an  observational  experiment  was  made  aboard 
the  Naval  Postgraduate  School  research  vessel,  R/V  ACANIA  shown  in 
Figure  1.   This  experiment  was  conducted  in  Monterey  Bay  approximately 
two  miles  offshore.   The  experimental  location  provided  the  near-ocean 
environment  for  turbulence  data  collection  as  well  as  the  platform  for 
ship  to  shore  optical  propagation  data  collection.  Meteorological 
measurements  of  wind  and  temperature  were  conducted  at  three  levels. 

The  scope  of  this  study  concerned  the  examination  of  the  turbulent 
fluctuations  of  velocity  and  temperature  in  the  marine  environment. 
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Figure  1.   The  R/V  ACANIA 


12 


The  examination  included  an  evaluation  of  the  height  variation  in  the 

constant  flux  layer  of  (1)  £   (the  viscous  dissipation  of  turbulent 

2 
kinetic  energy)  and  (2)  C    (the  temperature  structure  parameter). 
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II.   THEORETICAL  CONSIDERATIONS 


A.   THE  TEMPERATURE-STRUCTURE  PARAMETER 

An  important  parameter  in  describing  optical  propagation  is  the  re- 

2 
fractive-index-structure  parameter,  C   .   According  to  Wyngaard,  et  al 

2 
(1971),  C    can  be  related  to  its  counter  part,  the  tempera ture-struc- 

2 
ture  parameter,  C   through  a  semi-empirical  theory.   The  temperature- 
structure  parameter  has  the  form 

CT2  =  (79  x  10"6  P/T2)"2  CN2  (1) 

2 

C    can  be  estimated  by  several  methods  of  measurement.   Two  of  these 

methods  will  be  discussed  here. 

At  a  separation  distance,  r  ,  of  the  order  of  inertial  subrange 
scales,  the  temperature  structure  function  in  a  locally  isotropic  field 
has  the  form 


[T(x)  -  T(x  +  r)]2  =  CT2  r2/3  (2) 

2 

C    can  be  found  by  making  simultaneous  measurements  of  temperature, 

at  two  positions  separated  by  the  distance  r  .   A  second  method  for 

2 
determination  of  C    is  by  its  relation  to  the  one  dimensional  tem- 
perature spectrum  which  has  the  inertial  subrange  form 

<(>T(K)  =  6  CT2  K"5/3  (3) 


where  3  is  a  constant  and  K  is  the  streamwise  component  of  wave 
number . 
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2 
Indirect  means  of  expression  C    in  terms  of  other  meteorological 

parameters  is  through  turbulence  similarity  theory.  Surface-layer- 
similarity  theory  predicts  that  when  the  statistics  of  the  mean  and 
turbulent  flow  fields  are  properly  non-dimensionalized,  they  become 
universal  functions  of  only  a  stability  parameter.  Wyngaard,  et  al 
(1971)  chose  to  use  the  Richardson  number,  Ri,  as  the  stability  para- 
meter in  doing  calculations  from  overland  data.  The  Richardson 
number  is  defined  by 


Ri  =  S/T  <d9/f  >  (4) 

(duVdz) 


where  g  =  acceleration  of  gravity,   T  =  mean  temperature,   d0/dz  = 
potential  temperature  gradient,  and  du/dz  =  mean  wind  speed  gradient. 
It  was  shown  by  Wyngaard  that  the  temperature  structure  function  could 
be  expressed  as  a  function  of  the  Richardson  number,  potential  tempera- 
ture gradient  and  height. 

CT2  =  Z4/3  (d6/dz)2  f  (Ri)  (5) 

The  results  of  Wyngaard 's  calculations  are  shown  in  Figure  2,  verifying 

the  similarity  theory  prediction  for  unobstructed  overland  data. 

2 
In  order  to  fully  understand  the  relation  of  C    to  turbulence, 

one  must  consider   e   and   £Q  ,  the  viscous  dissipation  of  turbulent 

kinetic  energy  and  the  viscous  dissipation  of  temperature,  through 

theoretical  relations.   Both  the  balance  expression  for  turbulent 

kinetic  energy  and  Kolmogorov's  hypotheses  must  be  examined.   These 

theoretical  considerations  are  reviewed  in  the  following  paragraphs. 
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Figure  2. 
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B.   BALANCE  EXPRESSION  FOR  TURBULENT  KINETIC  ENERGY 

The  turbulent  kinetic  energy  balance  expression  has  the  following 
general  form 


2 


~  V.      8u.      ^  v.v. 

|r(4-)  =  -V^T^  +  ^V^  &.     -e+-|-  tv.(-^-^)] 

dt   2        XT  dx.    T    X    i_        dx.  1    p     2 
J   J    o       3         J 


s       v.v. 


Neglecting  horizontal  eddy  transports,  energy  flux  in  the  vertical  can 
be  expressed  as  follows 


2 

V.  du. 


£(-f )  -  -v3v.    ^    +  A  v.9  5   -  e  -  4    [v3[(^)   +  f]]    (7) 

3     o       3        3  L  J 


The  four  terms  on  the  right  hand  side  of  the  equation  represent  pro- 
cesses which  will  cause  a  time  variation  of  turbulent  kinetic  energy. 
These  processes  are 

i)  Mechanical  production  of  turbulent  kinetic  energy  from  the 
mean  flow 

ii)   Archimedes  work  term  or  the  energy  released  due  to  the  effects 
of  buoyancy 

iii)   Viscous  dissipation  of  turbulent  kinetic  energy 

iv)   Pressure  work  and  flux  divergence  of  turbulent  kinetic  energy 
due  to  transfer  by  turbulent  motion 

Further  simplification  of  the  kinetic  energy  balance  expression  can  be 

made  by  making  several  definitions  and  assumptions: 

i)  Assuming  that  the  turbulent  kinetic  energy  is  in  a  statisti- 
cally steady  state,      

2 

It  4"'  -  • 
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ii)   Assuming  horizontal  homogeneity  (i.e.,  the  spatial  variations 
of  mean  fluctuations  are  negligible), 

i-.i.-  o 

3x   dy 
iii)   Specifying  the  co-ordinate  system  such  that 


u2  =  u3  =  0 


iv)  Assuming  the  divergence  terms  are  negligible 


h3  [*3  ^  +  f  i  - » 


The  turbulent  kinetic  energy  balance  expression  can  then  be  expressed 
in  the  following  simplified  form,  for  the  case  when  the  mechanical  pro- 
duction of  turbulent  kinetic  is  equal  to  the  viscous  dissipation  plus 
energy  released  or  gained  due  to  buoyance  effects 

35, 


i  =  £  -  #-  vo0  (8) 


—V  V 

13  9x0       T   3 
3       o 


In  cases  where  the  atmosphere  is  in  near-neutral  stability  conditions, 
the  second  term  on  the  right  hand  side  of  the  equation  is  negligible 
and  Equation  (8)  can  be  written  as 

U*   9x3  "  G  (9) 

2 


where  u^  =  -v..v_  .   The  following  expression  for  the  vertical  gradient 

of  the  mean  wind  also  applies  to  the  neutral  condition 

9u    u 

(10) 


9x    Kz 
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where  K  is  a  universal  non-dimensional  constant  (i.e.,  the  Von  Karman 
constant)  and  Z  is  height.   Combining  Equations  (9)  and  (10)  yields 
an  expression  to  estimate  uA  based  on  measurements  of  e  at  a  known 

height  Z 

3 
u* 

e=  -B  (11) 


C.   VARIATION  OF   E  WITH  HEIGHT 

Equation  (11)  can  be  used  to  verify  the  constant  flux  assumption. 
A  straight  line  plot,  with  a  -1  slope,  of  lne  against  InZ  implies 
u.  is  a  constant  with  height.  For  non-neutral  stratification  in  the 
atmosphere,  one  must  retain  the  buoyancy  term  in  Equation  (8) .  Equa- 
tion (11)  can  be  modified  as  follows  to  include  the  buoyancy  effects. 

3 

e  ■  w-  [h  (b  -  b  (12) 


or 


3 


u* 


E  =  KZ~  H2    (R±)]  (13) 


_  3    

where  L   =  -Tu^  /gK  v  8  (Monin-Obuhkov  Length) 


V.0  =  heat  flux 

Ri  =  Richardson  number 

(J»1  a  empirical  function  of  Z/L 

(j)„  ■  empirical  function  of  Ri 

Equation  (13)  can  be  written  in  the  following  form 

3 
u* 
In  e  =  -In  Z  +  In  <{>0(Ri)  +  In  (-— -)  (14) 
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From  overland  experiments,  various  expressions  for  <f>   and  <}>   have 
been  derived.   For  cases  of  neutral  stability,  very  little  variation 
from  a  -1  slope  is  noted  in  a  plot  of   In  e  versus  In  Z  .   In  cases 
of  instability,  an  effect  on  the  plot  of   In  e  versus  In  Z  is  to 
increase  the  slope  (in  a  negative  sense),  to  greater  than  -1.   Figure  3 
illustrates  these  effects. 

D.   VARIATION  OF   C  2  WITH  HEIGHT 

2 
Equation  (5)  expressed  C    as  a  function  of  the  Richardson  number, 

potential  temperature  gradient  and  height.   The  potential  temperature 

gradient  can  be  expressed  in  the  following  manner 


where        T^  =  v  6/Ku^ 


L   ■  Monin-Obuhkov  Length 

g.   =  empirical  function  of  stability 


The  Richardson  number  can  also  be  expressed  as  a  function  of  Z/L  in 
Equation  (5) 

Ri  =  g2  (f)  (16) 

Combining  Equation  (5)  and  Equation  (16)  results  in  the  following  ex- 
2 


pression  for  C 


CT2  =  T*2  Z~2/3  g3  (f)  (17) 
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where  g_   is  an  empirical  function  of  stability  derived  from  overland 
data.   For  cases  of  neutral  stability,   g,  =  1  ,  and 

2     2  -2/3 
cT*   -  V  Z  <18> 


2 
A  plot  of  In  C    versus  In  Z  would  yield  a  straight  line  with  a  -2/3 

slope. 

2 
Under  stable  conditions  Wyngaard  noted  that  C    decreased  with 

-2/3 
height  more  slowly  than  Z     .   For  unstable  conditions,  overland 

results  indicate  that  g„  approximates  the  function 


83  "  I  (L)_2/3  (19) 


Also,  for  unstable  conditions,  Equation  (17)  can  be  expressed  as 

Ci^|T^(.L)2/32-4/3  (M) 


2 
A  plot  of  In  C    versus  In  Z  would  yield  a  straight  line  plot  with 

a  slope  of  approximately  -4/3,  assuming  overland  results  are  valid 

2 
over  water.   Figure  4  shows  the  predicted  variation  of  C    with 

height  for  various  stability  conditions. 


E.   KOLMOGOROV * S  SPECTRAL  SIMILARITY  EXPRESSIONS 

2 
The  parameters  e  and  C    can  be  estimated  on  the  basis  of 

Kolmogorov's  expressions  for  universal  spectra  (in  the  inertial  sub- 
range) of  both  fluctuating  velocity  and  temperature  spectra.   The 
inertial  subrange  in  the  velocity  spectrum  represents  a  range  of  eddy 
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Figure  A.  An  Illustration  of  the  Effect  of  Stability  on  the  Variation 
of   C    with  Height 
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sizes  where  energy  is  transferred  from  the  production  region  to  the 
dissipation  region,  (i.e.,  from  lower  to  higher  wave  numbers). 
Figure  5  from  Lumley  and  Panofsky  (1964)  shows  the  separate  regions 
of  a  typical  turbulence  spectrum. 

According  to  Kolmogorov's  second  hypothesis,  the  following  func- 
tional relations  hold 

<|>(K)  =  C.  £2/3  K~5/3  (21) 

4>T(K)=  C2  £_1/3£e  K"5/3  (22) 

where  <|>(K)   and  <j>_(K)   represent  one  dimensional  velocity  and  tem- 
perature spectra  and  K  represents  the  wave  number.   C.   and  C„  are 

2  -1/3 

empirical  constants.   C    is  defined  to  be  e  £fl  ;  therefore, 

Equation  (22)  is  identical  to  Equation  (3) .   Equations  (21)  and  (22) 

2 
can  be  used  to  find  estimates  of  e     and  C    by  examining  respective 

velocity  and  temperature  spectra. 

Turbulent  measurements  are  made  at  a  fixed  point  in  the  flow  and 

hence  spectra  are  obtained  for  temporal  frequency   f  .   In  order  to 

use  Equations  (21)  and  (22),  time  and  space  scales  must  be  related 

using  Taylor's  (1938)  "frozen  turbulence"  hypothesis 

K  =  2iTf/u  (23) 

where  u   is  the  mean  wind  speed.   The  hypothesis  implies  that  turbu- 
lence remains  unchanged  during  the  time  required  for  it  to  sweep 
past  the  probe.   The  following  forms  of  Equations  (21)  and  (22)  feature 
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Figure  5.   A  Schematic  Drawing  Showing  Spectral  Transfer  and 

Illustrating  (A)  Energy  Containing  Region  (B)  Inertial 
Subrange  (C)  Dissipation  Range  and  (D)  Universal 
Equilibrium  Range 
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the  transformation  from  frequency  to  wave  number 

.2/3  -2/3 


f  «J»(f)  -  K  (j»(K)  =  C±   e  ,J  K  ,J  (24) 

f  <j)T(f)=  K  <|)T(K)-  C2  CT2  K~2/3  (25) 


A  second  method  for  estimating  e  results  from  examining  the 
upper  limit  of  the  inertial  subrange  of  a  velocity  spectrum.   In 
statistical  equilibrium,  a  spectrum  of  turbulence  is  independent  of 
the  source  of  the  turbulence  and  depends  only  on  £  ,  y     and  K 
where  y     is  the  kinematic  viscosity.   The  following  expression  is 
applicable  for  the  inner  scale 

lo  =  (Y3/e)1/4  (26) 

where   (lo)    is  proportional  to  the  dissipation  wave  number.   The 
dissipation  wave  number  can  be  obtained  from  spectra  which  show  scale 
separation  between  the  inertial  subrange  and  viscous  dissipation 
regions,   e  can  then  be  determined  from  Equation  (26). 
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III.   OBSERVATIONS 

A.  DATA  COLLECTION 

Observational  experiments  were  made  from  the  R/V  ACANIA  anchored 
in  Monterey  Bay  approximately  two  miles  offshore.   Observations  were 
made  at  three  levels  on  the  forward  mast  and  bow  of  the  ACANIA. 
Figure  6  shows  the  mounting  arrangements  during  the  experiments. 

B.  EQUIPMENT 

1.  Fluctuating  Velocity  and  Temperature  Measurements 

Velocity  sensors  were  (TSI)  Model  1210  "hot  wire"  probes  with 
sliding  support  shields.   Calibration  of  the  sensors  was  made  prior 
to  the  experiments.   Band  filtering  was  applied  to  the  signal  with  cut- 
off frequencies  of  2  hz  and  2  khz  respectively.   The  signals  were  then 
amplified  and  recorded  on  a  14  channel  Sangamo  Model  3562  FM  recorder. 

Temperature  sensors  used  were  (TSI)  Model  1210  platinum  wire 
probes  with  sliding  support  shields.   Electronic  conditioning  similar 
to  that  applied  to  the  velocity  signals  was  applied  to  the  temperature 
signals.   The  temperature  data  was  analog  recorded  simultaneously  with 
the  velocity  data. 

2.  Mean  Velocity,  Temperature  and  Humidity  Measurements 

Mean  velocity  measurements  were  made  at  each  level  using  C.  W. 
Thornthwaite  cup  anemometers.   The  cup  revolutions  were  counted  by 
Hewlett  Packard  5221A  electronic  counters.   Direct  readout  of  mean  wind 
data  were  then  made  and  recorded  in  a  log.   Mean  temperature  measure- 
ments at  each  level  were  made  using  Hewlett  Packard  quartz  thermometers. 
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Figure   6.      Mounting  Arrangements 
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The  signals  were  passed  to  a  printer,  allowing  the  temperature  data  to 
be  printed  out  in  a  time  series.   Humidity  data  were  collected  in  a 
fashion  analogous  to  temperature  data  using  a  Hydrodynamics  Digital  II 
system  and  a  printer. 

Figures  7  and  8  show  the  electronics  involved  in  processing 
both  mean  and  fluctuating  data  in  the  ACANIA's  laboratory. 

C.   CALIBRATION 

1.   Calibration  of  Hot  Wire  Sensors 

Velocity  fluctuation  measurements  were  made  with  hot  wire  sensors. 
These  velocity  sensors  were  calibrated  before  each  experiment  using  a 
manometer  and  TSI  1125  Calibrator  in  accordance  with  procedures  outlined 
by  Johnston  (1974).   For  any  given  mean  wind  speed,  there  is  a  corres- 
ponding voltage  output.   The  relationship  between  wind  speed  and  voltage 
is  described  by  the  following  equation 

v2  =  a(u)1/2  +  b  (27) 

where   v   is  voltage  and  u   is  the  wind  speed.   Constants   a  and  b 

are  calibration  curve  slope  and  intercept  respectively,  and  are  determined 

experimentally . 

To  analyze  the  recorded  velocity  data  it  was  necessary  to  deter- 
mine a  calibration  factor  such  that 

(28) 
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Figure  7.   Equipment  Used  Aboard  R/V  ACANIA  to  Obtain  Fluctuating 
Velocity  and  Temperature  Data. 
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ts 

Figure  8.   Equipment  Used  Aboard  R/V  ACANIA  to  Obtain  Mean  Velocity 
and  Temperature  Data. 
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where  c  =  calibration  factor,   v'  ■  voltage  fluctuation  and  u'  = 
velocity  fluctuation.   Differentiation  of  Equation  (27)  yields 

[±JH ]  v«  =  u.  (29) 

3 

where  c  ,  the  calibration  is  then  determined  by  the  relation 

c  =  [^L—3  (30) 

For  analysis,  determination  of  the  calibration  factor  consisted  of 
(1)  obtaining  the  slope  of  the  calibration  curve  (2)  selecting  a  mean 
wind  speed  value  observed  during  the  experiment  and  (3)  entering  the 
calibration  curve  with  this  mean  wind  value  to  determine  a  correspond- 
ing voltage.   Figure  9  is  an  example  of  an  actual  hot  wire  calibration 
curve. 

2.   Calibration  of  Platinum  Wire  Sensors 

Temperature  fluctuation  measurements  were  made  using  platinum 
wire  sensors.   A  small  temperature  fluctuation  causes  a  corresponding 
resistance  change  in  the  platinum  wire.   The  resistance  of  a  given  wire 
can  be  related  to  temperature  through  the  following  equation 

R  =  Ro  (1  +  a  (T-To))  (31) 

where  R  =  resistance,   T  =  temperature,   To  =  0°C,   Ro  =  resistance 
at  0°C,  and  a  =  a  constant  for  platinum.   The  resistance,  in  ohms,  can 
be  measured  at  any  given  temperature  and  thus  Ro   determined  for  each 
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Figure  9.   Hot  Wire  Calibration  Plot 
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wire.   Differentiating  Equation  (31)  results  in 

dR  =  a  Ro  dT  (32) 

or 

R*  =  a  Ro  T'  (33) 

which  relates  a  small  temperature  change  to  a  small  change  in  the 
resistance  of  the  wire.   During  analysis,  resistance  fluctuations  are 
recorded  as  voltage  fluctuations,  through  the  electronics  involved. 
A  transfer  function,  which  relates  resistance  to  voltage,  must  be 
determined  during  calibration,  i.e., 

v'  =  HR'  (34) 

H  ,  the  transfer  function  can  be  found  by  substituting  a  variable 
resistor  in  place  of  the  wire  and  recording  voltage  changes  as  a 
function  of  corresponding  resistance  changes.   By  plotting  these  points 
on  a  graph,  the  resultant  slope  will  determine  the  transfer  function 
H  .   Figure  10  is  an  example  of  an  actual  calibration  curve  used  during 
one  of  the  experiments.   Substituting  Equation  (34)  into  Equation  (33) 
results  in  an  expression  to  relate  recorded  voltage  fluctuations  to 
temperature  fluctuations 

T'  =  <idb> v'  (35) 


or 


T'  =  C„v'  (36) 

n 


where  C   is  the  calibration  factor. 
H 
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Figure  10.   Platinum  Wire  Transfer  Function  Plot 
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IV.   DATA  REDUCTION  AND  ANALYSIS 

A.  PRELIMINARY  ANALYSIS 

Data  period  selection  was  based  on  evaluations  of  strip  chart 
records  of  the  signals.   Criteria  for  selection  of  data  periods  con- 
sisted of  (1)  signal  strength  (2)  possible  averaging  time  and  (3) 
availability  of  mean  profiles.   Typical  signals  meeting  the  above 
criteria  are  shown  in  Figure  11. 

B.  DIGITAL  METHODS 

1.  Analog  to  Digital  Conversion 

Data  from  selected  periods  were  prepared  for  further  analysis 
by  an  analog  to  digital  conversion  of  the  magnetic  tapes  as  described 
by  Johnston  (1974).   The  analog  to  digital  process  was  performed  on 
the  hybrid  computer  system  located  on  the  fifth  floor  of  Spanagel  Hall 
at  the  Naval  Postgraduate  School.   The  hybrid  system  consists  of  an 
analog  computer,  COMCOR  Ci  5000,  interfaced  electrically  with  a 
digital  computer,  XDS  9300. 

2.  Seven  Track  to  Nine  Track  Conversion 

The  analog  to  digital  process  yields  seven-track  octal  based 
records  written  on  a  digital  tape.   Analysis  of  the  data  was  then  per- 
formed on  an  IBM-360/67  computer.   This  computer  requires  nine-track 
hexadecimal  based  digital  records.   Conversions  from  seven-track  to 
nine-track  tapes  were  made  with  procedures  described  in  Naval  Post- 
graduate School  Computer  Center  Technical  Note  #0211-08.   Two  programs 
used  for  this  conversion  were  TAPEOUT  and  CONVERT. 
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3.   Digital  Analysis  Programs 

Spectral  analysis  of  the  digital  tapes  was  made  using  FTOR 
(Fast  Fourier  Transform)  and  SCOR  (Spectral  Analysis) ,  programs  which 
are  available  in  the  IBM  360/67  library.   Descriptions  of  the  FFT 
program  package  have  been  described  by  McKendrick  (1972)  and  Johnston 
(1974).   The  output  of  th4  SCOR  program  consists  of  a  tabular  computer 
printout  as  shown  in  Figure  12  and  a  graphical  plot  of  the  spectrum 
versus  frequency  as  shown  in  Figure  13. 

C.   SCALING  PROCEDURES 

During  observational  and  analysis  procedures,  several  gains  and 
scaling  factors  were  applied  to  the  recorded  signals.   Prior  to  record- 
ing the  data,  an  amplitude  gain  was  applied  to  the  signals.  During 
the  analog  to  digital  process  an  additional  gain  was  applied  to  the 
signal  along  with  a  scaling  factor.   The  scaling  factor  resulted  from 

the  consideration  that  the  maximum  voltage  input  to  the  COMCOR  Ci  5000 

23 
is  100  volts,  and  this  corresponds  to  a  digital  maximum  of  2   on  the 

XDS  9300.   Applying  the  respective  gain  and  scaling  factors  resulted 

in  a  signal  input  to  the  IBM  360/67  as  follows 

v(volts)  x  Gl  x  G2  x  SC  =  v(volts  IBM  360)        (37) 

where  v  =  signal  from  velocity  or  temperature  unit 

Gl  =  amplification  prior  to  recording 

G2  =  amplification  during  A-D  process 

SC  =  analog  to  digital  scaling  factor 
The  spectral  analysis  programs  yielded  plots  of  the  energy  density 
function  4>(f)  >  versus  frequency   f  .   The  corresponding  units 
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Figure  13.   SCOR  Spectral  Plot 
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2 
associated  with  cf>(f)   and  f  were  VOLT  SEC  and  1/SEC  respectively. 

Multiplying  each  point  of  the  energy  density  function  by  its  correspond- 

2 
ing  frequency  <}>(f)f  results  in  a  quantity  with  the  units  of  VOLTS 

Taking  the  square  root  of  this  quantity,  one  can  relate  VOLTS   to 

VOLTS'  (IBM  360)   through  Equation  (37).   Final  scaling  was  accomplished 

by  using  Equations  (29)  and  (35)  which  relates  volts  to  either  velocity 

or  temperature  fluctuations  respectively. 
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V.   INTERPRETATION  OF  SPECTRA 


A.   SINGLE  POINT  METHOD 

The  Single  Point  Method  is  perhaps  the  most  straightforward  method 

2 
used  to  obtain  e     and  C„   from  velocity  and  temperature  spectra. 

From  the  spectrum,  a  wave  number,   K  ,  within  the  inertial  subrange  is 

2 
chosen.   Using  Equations  (21)  and  (22),  calculations  of  e  and  C 

can  then  be  made  directly. 


B.   INTERCEPT  METHOD 

A  second  method  described  by  Johnston  (1974)  can  also  be  used  to 

2 
evaluate  e  and  C    from  velocity  and  temperature  spectra.   This 

method  has  an  advantage  over  the  single  point  method  in  that  many  points 
are  used  to  evaluate  the  spectral  estimate.   This  is  accomplished  by 
taking  the  average  slope  through  many  points  and  determining  an  inter- 
cept value  as  discussed  in  the  following  paragraphs. 
1.   Determination  of   £ 

Equation  (24)  can  be  written  in  the  following  form 

In  K<f>(K)  =  -2/3  In  K  +  In  C^2'3  (38) 


Plotting  In  K$(K)   against  -2/3  In  K  results  in  a  curve  with  a  slope 

2/3 
of  +1  and  an  intercept  of   In  C,£     .   At  the  point  of  intercept, 

In  K  equals  zero  and 

In  K<J>(K)  =  In  C^c2'3  (39) 
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Letting  A  equal  In  K<j)(K)  ,  the  intercept  value,   £   then  becomes 

e  =  (|V/2  (40) 


By  determining  the  intercept  value,   £  can  be  estimated.   Figure  14 

shows  a  typical  plot  of  this  type  using  velocity  data  at  three  levels. 

2 
2.   Determination  of  C 

2 
A  similar  procedure  can  be  used  to  determine  C    .   Equation 

(25)  can  be  written  as  follows 

In  K*  (K)  =  -2/3  In  K  +  In  C2  C^,2  (41) 


Using  the  same  procedure  as  described  in  the  determination  of   e  , 

2 
C„   can  be  evaluated  at  the  point  of  intercept  in  the  following  manner 

where  B  =  In  K<}>  (K)   at  the  ordinate  intercept.   Figure  15  is  an 
example  of  an  intercept  plot  using  temperature  data  at  two  levels. 


C.   INNER  SCALE  METHOD 

A  third  method  of  estimating   £  was  to  consider  the  inner  scale. 
Equation  (26)  describes  an  expression  relating  £  to  the  separation 
between  the  inertial  subrange  and  dissipation  region  of  the  spectra. 
By  choosing  the  wave  number  at  which  the  dissipation  region  occurs  in 
the  velocity  spectra,   £   can  be  evaluated.   This  method  was  quite  sub- 
jective with  regard  to  determining  the  exact  point  of  scale  separation. 
However,  the  results  were  qualitatively  good  when  compared  to  the 
intercept  method. 
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Figure  14.   Intercept  Plot  for  Velocity  Data  at  Three  Levels, 
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Figure  15.   Intercept  Plot  for  Temperature  Data  at  Two  Levels. 
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VI.   RESULTS 

A.  SUMMARY  OF  PERIODS  EXAMINED 

Fluctuating  velocity  and  temperature  data  were  collected  from  three 
levels  aboard  the  R/V  ACANIA  on  18  January  1974.   Selected  periods  for 
analysis  were  chosen  for  three  time  intervals  between  1834  and  1932. 
Table  I  lists  periods  examined,  heights  at  which  data  were  collected, 
parameters  measured  and  stability  conditions  for  each  period. 

B.  VELOCITY  RESULTS 

Velocity  spectra  were  obtained  for  simultaneous  measurements  at  three 
levels  during  each  period.   From  these  spectra  e     and  uA  were  computed 
by  methods  described  in  Sections  V-B-l  and  V-C.   Table  II  lists   £   and 
u.   calculations  for  each  of  these  periods. 

The  original  experimental  plans  were  to  measure  mean  wind  data  at 
each  of  the  levels.   Due  to  problems  encountered  with  the  read-out  counter, 
the  wind  data  appeared  to  be  in  error.   The  problem  was  resolved  by  use  of 
mean  wind  measurements  at  one  level,  the  ACANIA 's  aerovane,  and  then  esti- 
mating a  vertical  wind  field  based  on  the  Monin-Obuhkov  prediction.   A 
Fortran  computer  program,  PROPAR,  described  by  Cardone  (1969)  was  used 
to  estimate  the  vertical  wind  field.   This  program  uses  inputs  of  (1) 
mean  wind  speed  (2)  anemometer  height  (3)  air-sea  temperature  difference 
and  (4)  thermometer  height  to  compute  the  output  parameters.   These  com- 
puted parameters  included  (1)  u^.  =  friction  velocity  (2)  Zo  =  roughness 
parameter  and  (3)  L  =  Monin-Obuhkov  Length.   Mean  wind  was  calculated 
at  each  level  using  the  following  equation 
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TABLE  I 

Data  Periods  Considered  in  This  Study 

18  January  1974 


Time  Height  Parameter  Stability 

(local)  (cm)  (Ri) 


1834-1854  394.0  U,T,Q,U',T'  -0.136 

(level  1) 

658.0  U,T,Q,U',AT' 

(level  2) 

822.0  UjfjQ.U'.T1 

(level  3) 


1855-1913  394.0  U,f,Q,U',T'  -0.293 

(level  1) 

658.0  U,T,Q,U',AT* 

(level  2) 

822.0  U,T,Q,U*,T' 

(level  3) 


1914-1932  394.0  U,T,Q,U*,T'  -0.923 

(level  1) 

658.0  U,T,Q,U',AT' 

(level  2) 

822.0  U,T,Q,U',T' 

(level  3) 
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TABLE  II 

e  and  Momentum  Flux  Calculations  for  18  January  1974  using 
the  Zero  Intercept  Method 


Height                          e  u^ 

(cm  sec  )  .      -1. 

(cm  sec  ) 

Time  1834-1854 

394.0                        83.43  23.60 

658.0                        54.74  24.33 

822.0                         22.83  19.57 

Time  1855-1913 

394.0                       342.96  37.80 

658.0                         79.21  27.51 

822.0                         36.72  22.93 

Time  1914-1932 

394.0                        112.28  26.05 

658.0                         13.67  15.32 

822.0                         22.22  19.40 
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u  =  ir  (ln  <i>  -  *i  <fb)  <43> 


where  i^1   is  an  empirical  function  of  Z/L  described  by  Paulson  (1970). 
The  purpose  in  defining  the  wind  field  distribution  was  to  obtain  advec- 
tion  speeds  for  relating  temporal  and  spatial  frequencies.   The  need  for 
wave  number  (spatial  frequency)  description  was  identified  in  Section 
II-E. 

Figures  16,  17  and  18  show  typical  velocity  spectra  at  each  level  for 
one  of  the  periods.   In  each  case  a  -5/3  slope  is  observed  for  the 
inertial  subrange  region.   At  higher  frequencies,  the  spectra  exhibit 
greater  slopes.   The  regions  of  different  slopes  represent  different 
scales  in  the  turbulent  regime.   The  scale  represented  by  greater  slope 
at  higher  frequencies,  is  the  dissipation  region.   Closer  investigation 
reveals  a  decrease  in  the  frequency  where  the  slope  changes  as  height 
increases.   However,  as  mentioned  earlier  and  as  can  be  seen  here,  the 
exact  point  of  scale  separation  between  the  inertial  subrange  and  dissi- 
pation region  is  difficult  to  determine. 

Results  on  e  versus  height  for  the  velocity  data  are  shown  in 
Figures  19,  20  and  21  where  In  e     is  plotted  against   In  Z  for  each 
data  period.   Negative  slopes,  greater  than  -1,  are  noted  in  each  case. 
As  stated  in  Section  II-C,  this  agrees  with  present  boundary  layer  theory 
for  unstable  atmospheric  conditions. 

The  deviation  of  the  slope  from  -1  was  examined  to  see  if  it  could 
be  accounted  for  by  the  influence  of  the  unstable  stratification.   The 

plots  were  corrected  for  stability  using  Equation  (13)  where  <}>„  = 

-1/4 
(1  -  12Ri)      .   A  slight  shift  toward  a  -1  slope  was  observed;  however, 
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Figure  16.   Velocity  Spectrum  Level  1  for  1855-1913,  18  Jan  74. 
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Figure  17.   Velocity  Spectrum  Level  2  for  1855-1913,  18  Jan  74. 
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Figure  18.   Velocity  Spectrum  Level  3  for  1855-1913,  18  Jan  74, 
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Figure  19.      Measured  Values  of      e     vs  Height   for   1834-1854, 
18  Jan   74. 
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Figure  20.   Measured  Values  of   e  vs.  Height  for  1855-1913,  18  Jan  74, 
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Figure  21.  Measured  Values  of  e  vs.  Height  for  1914-1932,  18  Jan  74. 
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complete  adjustment  to  a  -1  slope  was  not  evident.   Similar  results 
were  found  by  Johnston  (1974).  He  states  "The  fact  that  the  stability 
corrections  did  not  completely  restore  the  -1  slope,  even  in  the 
stable  case,  leads  to  the  examination  of  the  possible  effects  of  wind- 
wave  coupling." 

The  above  results  imply  that  overland  prediction  for  the  influence 
of  stability  are  not  completely  valid  over  water  and  that  wind-wave 
coupling  or  wave  related  velocity  fluctuations  might  be  the  reason. 
There  is  still  not  enough  experimental  and  theoretical  knowledge  to 
describe  the  waves'  influence  upon  the  data. 

Several  consistent  features  can  be  observed  in  the  results.   During 
the  first  two  periods  there  appears  to  be  a  bending  or  "kink"  in  the 
plots  such  that  they  are  concave  downward.   In  the  third  data  period, 
the  "kink"  is  concave  upward.   However,  strip  chart  records  show  a 
deterioration  of  the  velocity  signal  at  level  2  toward  the  end  of  the 
period.   This  apparently  resulted  in  low  estimates  of  e  for  that 
level.   It  is  not  felt  that  this  "kink"  is  due  to  the  ship's  influence. 
Welsh  (1974)  found  that  ship  related  motion  had  little  effect  on  turbu- 
lence measurements  at  these  frequencies.   Johnston  (1974)  did  not  find 
this  characteristic  bending  in  similar  plots  of  In  £  versus  In  Z 
for  unstable  conditions.   In  these  results,  both  the  intercept  method 
and  inner  scale  method  exhibit  the  "kink".   Based  on  only  two  periods 
of  analysis,  conclusions  as  to  the  cause  of  this  feature  cannot  be 
asserted. 

C.   TEMPERATURE  RESULTS 

Fluctuating  temperature  data  were  also  obtained  at  three  levels  for 

2 
each  period.   At  levels  1  and  3,   C    was  determined  from  the 
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one-dimensional  temperature  spectrum  as  described  by  Equations  (3)  and 

(22).  At  the  mid  level,  two  probes  were  separated  by  r  =  10  cm  and 

2  2 

C    determined  by  using  Equation  (2).   Table  III  lists  the  C    calcu- 
lations for  each  of  the  periods.   Temperature  spectra  at  three  levels 
are  shown  for  one  of  the  periods  in  Figures  22,  23  and  24.   These  are 

typical  of  the  spectra  for  the  other  periods  examined. 

2 
Results  of  C    versus  height  are  shown  in  Figures  25,  26  and  27, 

2 
where  In  C    is  plotted  against  In  Z  .   The  general  slopes  of  these 

plots  agree  reasonably  well  with  the  -4/3  prediction  for  unstable  con- 
ditions discussed  in  Section  II-D.   Periods  1  and  2  both  show  close 
agreement  with  the  -4/3  prediction.   However,  period  3  shows  a  much 
1  greater  deviation  from  -4/3.   Comparing  Richardson  numbers  from  Table  I, 
it  is  noted  that  period  3  was  much  more  unstable  than  the  first  two 
periods.   One  might  expect  slightly  greater  negative  slopes  for  greater 


instability;  however,  -4/3  should  be  a  limiting  slope  if  overland  results 

2 
"T 


are  valid  over  water.   Since  Cm   is  a  function  of  both  e  and  e0  , 


the  possible  wind  wave  effects  on   £   discussed  in  the  previous  section, 

2 
would  also  influence  C    .   This  is  a  possible  explanation  for  any 

2 
deviation  in  C    from  a  -4/3  slope. 

2 
Closer  examination  of   In  C    versus   In  Z   reveals  "kinking"  or 

bending  similar  to  that  discussed  in  the  velocity  results  section.   In 

both  periods  1  and  2,  this  bending  phenomenon  is  concave  downward  as  it 

was  for   £   calculations.   During  period  3  the  general  appearance  of  the 

2 
kink  in  the   C    plot  approximates  the  appearance  of  the  kink  in  the   £ 

plot,  disregarding  any  influence  of  a  deterioration  in  the  velocity 

signal  during  that  period.   In  this  case  however,  the  bending  is  in 

the  opposite  direction. 
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TABLE  III 

2 

C    Calculations  for  18  Jan  74  using  the  Zero 

T   Intercept  Method 


Height 
(cm) 


(°c2 


-2/3, 
cm    ) 


394.0 
658.0 
822.0 


Time  1834-1854 


.00312 
.00239 
.00168 


394.0 
658.0 
822.0 


Time  1855-1913 


.00600 
,00558 
.00428 


394.0 
658.0 
822.0 


Time  1914-1932 


.00934 
.00416 
.00393 
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Figure  22.   Temperature  Spectrum  Level  1  for  1855-1913,  18  Jan  74. 
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Figure  23.  Differential  Temperature  Spectrum  for  Fluctuating 
Temperature  Data  Separated  by  r  =  10  cm,  Level  2 
for  1855-1913,  18  Jan  74. 
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Figure  24.   Temperature  Spectrum  Level  3  for  1855-1913,  18  Jan  74, 
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Figure  25.  Measured  Values  of  C    vs.  Height  for  1834-1854 
18  Jan  74. 
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Figure  26.   Measured  Values  of  C    vs.  Height  for  1855-1913, 
18  Jan  74. 
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Figure  27.   Measured  Values  of  CT   vs.  Height  for  1914-1932 
18  Jan  74. 
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The  above  results  give  further  support  to  the  consideration  that 
the  bending  is  actual  and  not  a  result  of  random  errors  within  the 
system,  since  these  spectra  were  obtained  from  two  independent  sources, 
velocity  and  temperature  measurements.  However,  no  further  explana- 
tion as  to  the  cause  of  the  bending  is  offered  due  to  the  limited 
amount  of  data  considered. 
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VII.   CONCLUSIONS  AND  RECOMMENDATIONS 

2 
Observed  height  variations  of  e  and  C    tend  to  support  present 

formulations  and  prediction  techniques.   However,  several  exceptions 

were  noted  and  will  be  summarized  in  the  following  paragraphs. 

For  unstable  atmospheric  conditions  slopes  of  In  e  versus  In  Z 

were  greater  than  -1.   When  these  slopes  were  corrected  for  stability, 

there  was  a  tendency  for  the  slopes  to  shift  towards  a  -1  slope.   The 

observation  was  made  that  complete  adjustment  to  a  -1  slope  was  not 

evident,  suggesting  other  factors  such  as  wind  wave  effects  may  be 

important  in  near-ocean  boundary  theory.   The  observation  was  also  made 

that  In  e     versus   In  Z  was  not  a  linear  plot.   There  appeared  to  be 

bending  or  "kinking"  which  is  not  predicted  in  present  formulations. 

No  definite  reason  was  given  to  account  for  the  apparent  kinking  due 

to  the  limited  amount  of  data  examined. 

2 
A  plot  of  In  C    versus  In  Z  resulted  in  -4/3  slopes  for  unstable 

conditions.   Some  slight  deviations  were  noted  and  attributed  to  other 

factors  mentioned  in  the  previous  paragraph.   The  characteristic  "kinking" 

2 

or  bending  in  the   In  C    versus   In  Z   plots  was  also  evident. 

Spectral  analysis  appears  to  be  a  satisfactory  method  for  determining 

2 
e  ,  momentum  flux  and  C    estimates  over  the  ocean.   However,  digital 

methods  used  are  very  time  consuming  due  to  the  amount  of  tape  handling 

and  large  computer  programs  used  in  spectral  analysis.   A  real-time 

spectrum  analyzer  would  increase  the  amount  of  data  that  one  could 

examine  due  to  the  time  saved  in  the  analysis.   It  is  suggested  that 
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subsequent  fluctuating  velocity  and  temperature  experiments  aboard  the 
ACANIA  be  examined  using  a  spectrum  analyzer  in  the  analysis  process. 
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